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POSTSEISMIC DEFORMATION DUE TO SUBCRUSTAL VISCOELASTIC 
RELAXATION FOLLOWING DIP SLIP EARTHQUAKES 


Steven C. Cohen 
ABSTRACT 

The defonnation of the earth foliowinj! u dip-slip earthquake is calculated usinji a 
three-iaver rheological model and finite element techniques. The three layers are an 
elastic upper lithosphere, a standard linear solid lower lithosphere, and a Maxwell 
viscoelastic asthenosphere-a model pawiously analyzed in the strike-slip case (Cohen. 
l*J8l. l‘J82). Attention is focused on tlie magnitude of the postseismic subsidena 
and the width of the subsidence zone that can develop due to the viscoelastic response 
♦o coseismic reverse slip. Detailed analysis for a fault extending from the surface to 
15 km with a 45° dip reveals that postseismic subsidence Is sensitive to the deptli to 
the asthenosphem but is only weakly dependent on lower lithosphere depth. Tlie greatest 
subsidence occurs when the elastic lithosphere is about 30 km thick and the asthenosphere 
lies just below this layer (asthenosphere depth = 2 times the fault depth). Tlie 
extremum in the subsidence pattern occurs at about 5 km from the surface trace of the 
fault and lies over tlie slip plane. In a typical case after a lime t = 30r (t = Maxwell 
time) following the earthquake the subsidence at this point is bO <• ot !lie coseismic 
uplift. Unlike the horizontal deformation following a strike-slip earthquake, significant 
vertical deformation due to a.sthenosphere fiow persists for many limes r and the 
magnitude of the vertical defonnation is not necessarily enhanced by having a partially 
relaxing lower lithosphere. Hie width of the postseismic subsidence zone also depends 
strongly on the depth to the asthenosphere and shows little sensitivity to the presence 
of a slowly relaxing lower iithosphere. The model fits postseismic rebound data from 
the I8db Rilai-u. Japan e irthquake. but the fit is not substantially better than that 
obtained with a simpler two layer model. Previous estimates of the values of 
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asthenosphere depth (~30-40 km) and viscosity poise) are confined by the 

three layer calculations. 
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POSTSEISMIC DEFORMATION DUE TO SUBCRUSTAL VISCOELASTIC 
RFT AXATION FOLLOWING DIP-SLIP EARTHQUAKES 


INTRODUCTION 

During the past few years many ps^pers have been published on models of crustal deformation 
during the earthquake cycle (e.g. Fitch and Scholz, 1971 ; Bischke, 1974; Rundle and Jackson, 
1977; Scholz and Kato, 1978; Savage and Prescott, 1978; Spence and Turcotte, 1979; Melosh and 
Fleitout, 1982). Among the mechanisms proposed to account for various attributes of the temp- 
oral-spatial deformation patterns observed after earthqua!:es are surface afterslip, aseismic slip at 
depth (Thatcher, 1974), viscoelastic flow of subsurface layers of the earth (Nur and Mavko, 1974; 
Lehner and Li, 1981 ;Matsu’ura, 1981 ;Cohen, 1981), and various other forms of surface and sub- 
surface yielding (Wahr and Wyss, 1980; Yang and Toksoz, 1981 ; Mavko, 1981 ). Comparisons of 
model calculations with geodetic measurements have revealed substantial agreement in a number 
of cases (e.g. Thatcher and Rundle, 1979; Brown et. al., 1977). Nevertheless, numerous questions 
remain unresolved about the most appropriate models to be used to represent the postseismic 
response to specific earthquakes. Many geodetic observations can be explained equally well by 
competing models and the qualitative and quantitative features of individual models sometimes 
show substantial variation depending on the choice of poorly constrained model parameters. 
Despite these problems progress has been made in putting numerical bounds on key descriptive 
parameters such as slip magnitude and depth, viscoelastic layer depths, and viscosity by perform- 
ing detailed calculations with several models With this in mind I have studied several layered 
viscoelastic models wherein layers of varying rheological constitutive laws and varying viscosity 
are distributed by depth (Cohen 1980, 1981, 1982). For example in this note and two 
previous papers an elastic upper lithosphere lies over a partially viscoelastic lower lithosphere, 
which in turn lies ever an even more fluid asthenosphere. This is similar but not identical to the 
situation where an elastic lithosphere lies over a layered viscoelastic asthenosphere. A somewhat 
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different layering can occur if the lithosphere itself can be divided into an elastic upper crust, a 
decoupling intracrustal low viscosity zone, and an elastic lower lithosphere (Turcotte et. al. 1981). 
The purpose of the present paper is to complete my analysis of multilayer viscoelastic models by 
extending earlier calculations of postseismic rebound from strike-slip earthquakes to the case of 
dip-slip events. 

MODEL DESCRIPTION AND RESULTS 

The model I wish to consider divides the earth into three layers. All the layers are elastic in 
their bulk behavior. The top layer (upper lithosphere) is elastic in shear. Lying within this layer is 
a dislocation plane (fault) that dips an angle Q from the horizontal and extends from an upper 
depth, d, to a lower depth, D. The layer thickness is Hj and 0 < d < D < Hj . The second layer 
(lower lithosphere) begins at depth Hj and extends to depth H 2 ; its thickness is H 2 - Hj . This 
layer behaves as a standard linear solid in shear. The third layer (asthenosphere) begins at depth 
H 2 . The layer acts as a Maxwell substance in shear. The rationale for this model is discussed in 
Cohen(1982). For the calculations reported here, the elastic moduli of the various layers are equal; 
the shear modulus has a value 3x10^^ dyne/cm^ and the bulk modulus the value 8.33 x 10^ ^ 
dyne/cm^. Viscosities of 17 ^ = 5 x 10*^ poise and tjj = 1 x 10^^ poise are assumed for the astheno- 
sphere and lower lithosphere respectively (Cohen, 1982). The asthenosphere extends to a depth of 
300 km. The other parameters that enter the model are the fault depths d and D, layer depths, 
and H 2 , and dip angle 6 . The results depend on the coseismic fault slip, the time elapsed since the 
earthquake, and the location. The general differential equations for the displacements have already 
been published (Cohen, 1982); here we present results using a plane strain calculation. The plane is 
a cross section through the earth; the fault is a dipping line segment m this plane. Slip is parallel to 
and along this line segment. The differential equations are solved by finite element techniques. 

In Figure 1 I show an example of the vertical and horizontal displacements as a function of 
distance from the surface trace of the fault and of time after the earthquake. Qualitatively many 
of the features are similar to those observed in elastic lithosphere-viscoelastic asthenosphere 
(2-layer) calculations. The dominant feature of the vertical displacement pattern is the broad 
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postseismic subsidence over the fault and uplift further away. This general shape for the vertical 
displacements is maintained over a wide range of layer depths and dip angles when the fault is short 
enough that it does not rupture nearly all the lithosphere (Melosh, 1982). It is convenient for 
presenting figures to normalize the results to the coseismic slip, Sq and to take the unit of logth 
as the fault depth, D (where for the res ilts I will show, in Figures 3 and 4, d = 0). The unit of time 
is the Maxwell time of the asthenosphere r = where ri^ and are the viscosity and rigidity 
respectively. Figure 2 shows an example of how the magnitude of the subsidence, A and the 
width of the subsidence zone, W, vary with time. As expected, A changes most rapidly immedi- 
ately after the earthquake when the asthenosphere begins to relax the high coseismic stress change. 
As stress relaj.ation proceeds, the subsidence rate decreases. The width of the subsidence zone 
generally grows at a decreasing rate with time. However, for a short time after the earthquake, the 
subside*''^e zone width decreases slightly before the longer term growth becomes dominant. The 
widi;h of *' v' subsidence zone is much broader than the coseismic uplift zone and in fact is larger 
than the depth to the asthenosphere, W>H 2 >Hj>D. It is interesting to focus attention on how 
the subsidence varies with layer depths. Figure 3 shows the value of the postseismic vertical dis- 
placement, A V, (t = 30 t) as a function of the lower lithosphere and asthenosphere depths. The dis- 
tance chosen, X = -1/3 D, is near the extremum in the subsidence pattern. The interesting feature 
of Figure 3 is that A V fails to behave monotonically for either an increase in H 2 or H j with the 
other held fixed. For H 2 fixed, A V increases as increases until H| > 2D then AV decreases. 
Similarly for Hj fixed AV increases with H 2 until H 2 > 2D then it decreases. Stated in another 
way: when the upper lithosphere thickness is more than 2D, subsidence is greatest when the lower 
lithosphere is vanishingly thin, when the upper lithosphere thickness is greater than 2D, subsidence 
is maximized when the lower lithosphere thickness, T = 2D - Hj . Similarly when the depth to 
tne asinenosphere is less than 2D, subsidence is greatest when the lower lithosphere vanishes, but 
when the asthenosphere depth is greater than 2D, subsidence is maximum for T = H 2 - 2D, The 
greatest subsidence occurs when Hj = H 2 = 2D. The condition for maxii.;um subsidence H 2 = 2D 

has been verified at other dip angles. The magnitude, AV = 0.26, is approximately 60% of 
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the coseismic uplift in the ^.resent case. Another interesting feature of Figure 3 is that 
AV changes by only modest amounts over a substantial range of Hj values. For example, when 
Hj = H 2 = 10/3 D, AV = 0.19, a vame not that much less than the extremum value AV = 0.26 at 
Hj = H 2 = 2D. Only wheri H 2 gets small and the postseismic patterns change from those for which 
subsidence is dominant to those having substantial uplift do rapid changes in AV arise. For com- 
parison with simpler two-layer models one can contrast the A V at various values of Hj with those 
where Hj = H 2 (i.e. viscoelastic lower lithosphere is vanishing thin). 

In Figure 4 I plot the width of the subsidence zone, W, as a function of Hj and H 2 (again at 
X = -1/3 D and t = 30 t). The most striking feature is the growth of W with increasing H 2 and the 
relative insensitivity to Hj . The width of the deformation zone is dominated by the depth to the 
fastest relaxing layer, the fact that a shallower zone is also relaxing (at a slower rate) has little impact 
on the spatial pattern. In the case of postseismic rebound to slip on a strike-slip fault, the deforma- 
tion due to asthenosphere flow is largely completed in a few Maxwell times. In the dip-slip case, 
however, it seems that the development of the deformation due to asthenosphere flow persists for 
much longer times. As a result, the existence of a relaxing viscoelastic rather than an elastic lower 
lithosphere does not necessarily enhance the long term subsidence. 

Figure 5 compares the temporal and spatial patterns of the vertical displacement as changes are 
made in the fault depth. Curves 1 and 2 refer to faults that rupture the surface but extend to dif- 
ferent depths; in tiiis case the postseismic displacements are very comparable. By contrast, curve 
3 refers to a fault that is buried and has the same length as that of curve 2 and the same lower depth 
as that of curve 1. The postseismic subsidence is much less in the buried fault case. Nevertheless, 
the shape of the postseismic deformation curves are similar in the three cases; certainly they are 
more similar than the coseismic curves. 

For a comparison of model predictions to observations I have chosen the 1896 Riku-u, Japan 
earthquake. As discussed in Thatcher, et. al. (1980) the subsidence data for this region is derived 
from surveys in 1900, 1934, and 1974. The model results, shown in Figure 6 are derived for the 
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case where both the lithosphere and asthenosphere have thicknesses of 20 km. They confirm 
Thatcher’s et.al. estimate of ~30 km for the lithosphere thickness and ~ 10^® poise for the 
viscosity in this region of Japan. Despite the long time span of the record, little information is 
derived about the lower lithosphere viscosity. Its thickness however, must not be substantially 
greater than 20 km. 

SUMMARY 

In the postseismic lebound to strike-slip earthquakes the effect of having a lower lithosphere 
which is viscoelastic rather than elastic is to increase both the horizontal displacement magnitude 
and the time duration of the deformation process. The present study has shown that the vertical 
displacements after dip-slip earthquakes have a more complicated dependence on the existence of 
a relaxing lower lithosphere. However, over a wide range of layer depths the surface displacements 
are far less influenced by lower lithosphere relaxation than by asthenosphere relaxation. Because 
the magnitude, duration, and width of the postseismic subsidence depend primarily on the depth 
and viscosity of the most fluid layer geodetic data can be used to deduce asthenosphere depth 
and viscosity, but little information is obtained about the lower lithosphere. 
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FIGURE CAPTIONS 


Figure I : Vertical displacements, V, and horizontal displacements, U, versus distance from fault 

at various times; coseismic displacements Vq, Uq at time tQ = 0'*’; postseismic displace- 
ments A V = V(tj) - Vq (tQ), AU = U(tj) - UQ(tQ) at tj = 1 • 10^ sec = 6 T, t 2 = 5 • 10^ 
sec = 30 T. Coseismic reverse slip = 1 m, dip angle = 45° ; fault depth = 0- 1 5 km 
asthenosphere viscosity = 5 • 10^^ poise; lower lithosphere viscosity 1 • 10^^ poise. 

Figure 2: Dependence of a subsidence extremum value, AV^.^, and subsidence zone width W on 

time, t, for a model with Hj/D = 4/3 ; H 2 /D = 8/3. Greatest subsidence occurs at 
X = -1/3 D. Sq is coseismic dip and r is Maxwell time of asthenosphere. Fault dip 
angle = 45°. 

Figure 3; Magnitude of postseismic subsidence, A V, as a function of the depth of the astheno- 
sphere, H 2 , and lower lithosphere, Hj . Distance from fault, X = -1/3 D, time after 
earthquake, t = 30 r. Other parameters are 0 = 45°, ~ 20. 

Figure 4: Width of postseismic subsidence zone, W, as a function of the depth of the astheno- 

sphere, H 2 , and lower lithosphere, Hj . Time after earthquake, t = 30 t. Other para- 
meters are 6 = 45“ . tjj/t?^ = 20. 

Figure 5; Vertical displacements for various fault depths; curve : d = 0 km, D = 15 km; curve 2; 
d = 0 km, D = 10 km; curve 3:d = 5 km, D = 15 km. 

Figure 6; Comparison of leveling data and model calculations for postseismic rebound to 1896 
Riku-u Japan earthquake a. 1934 1900 data; b. 1974- 1900 data. Survey results use 
a 3-point smoothing to data of Thatcher, et. a! ( 1980). Model parameters are coseismic 
slip = 2 m, fault depth = 0-15 km; dip angle = 45° ; asthenosphere depth = 40 km; 
asthenosphere viscosity = 5 x 10*^ poise, lower lithosphere depth = 20 km, lower 
lithosphere viscosity = 1 x lO"* poise. 
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Figure 1. Vertical displacements, V, and horizontal displacements, U, versus distance from fault at v^ous times; coseismic disr- 
placements Vq. Uq at time to = j*\ postseismic displacements AV = V(tj) - VgCto), AU = U(tj) ~ at tj = 1 • 10 sec = 

6 T, t2 = 5 • 109 sec = 30 t. Coseismic reverse slip = 1 m, dip i.igle = 45“ ; fault der^th = 0- i 5 km asthenosphere viscosity = 5 • 10* ^ 
poise; lower lithosphere viscosity 1 • 10^* poise. 
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Figure 2. P-pendence of a subsidence extremum value, AVj^, and subsidence zone width W on 
time, t, fora model with H,/D = 4/3;H2/D = 8/3. Greatest subsidence occurs at X = -1/3 D. S 
is "oseismic dip and t is Maxwell time of asthenospheie. Fault dip angle = 45 . 




Figure 3. Magnitude of postseismic subsidence, AV, as a function of the depth of the astheno- 
sphere, H2 , and lower lithosphere, Hj . Distance from fault, X = -1/3 D, time after earthquake, 
t = 30 T. Other parameters are d = 45°, tjj/tJj = 20. 
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FiKi."- 4 Width cf postseismic subsidence zone, W, as a function of the depth of the astheno- 
spLre. H2 . and lower lithosphere, Hj . Time after earthquake, t = 30 r. Other parameters are 

e = 20. 


11 







SURVEY DATA 
MODEL CALCULATIONS 


ORIGINAL PAGE fS 
OF POvOR QUALITY 



(WD) lN3lM30VTdSia 1V0lia3A 


X 



mo) iN3M30Vl«iSia 1V3I1U3A 


13 


Figure 6. Comparison of leveling data and model calculations for postseismic rebound to 1 896 Riku-u Japan earthquake a. 1934- 
1900 data; b. 1979-1900 data. Survey results use a 3-point smoothing to data of Thatcher, et.al. (1980). Model parameters are 
coseismic slip = 2 m, fault depth = 0-15 km; dip angle = 45° ; asthenosphere depth = 40 km; asthenosphere viscosity = 5 x lO^^, 
poise, lower lithosphere depht = 20 km, lower lithosphere viscosity = 1 x 1021 poise. 



